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Abstract 

The spectrum of D = 4 supersymmetric Yang- Mills quantum mechanics is com- 
puted with high accuracy in all channels of angular momentum and fermion number. 
Localized and non-localized states coexists in certain channels as a consequence of 
the supersymmetric interactions with flat valleys. All states fall into well identifiable 
supermultiplets providing an explicit realization of supersymmetry on the spectro- 
scopic level. An accidental degeneracy among some supermultiplets has been found. 
Regularized Witten index converges to a time-independent constant which agrees 
with earlier calculations. 
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1 Introduction 



In the present paper we report detailed studies of the supersymmetric Yang-Mills quantum 
mechanics (SYMQM) [1],[2]. The particular model addressed here results from the di- 
mensional reduction of the supersymmetric Yang-Mills field theory, with the SU(2) gauge 
group, from a four-dimcnsional space-time {D = 4) to a single point in space. It is a mem- 
ber of a family of quantum mechanical systems with the famous D = 10, SU(iV — > oo), 
SYMQM at its upper end. The latter, considered as a model of an M-theory [3] , attracted 
a lot of attention in recent years, see [4] and [5] for reviews and further references. For 
that reason we have launched a systematic, nonperturbative study of the whole family 
(varying D and A^) in an attempt to understand their global properties, and to develop 
adequate techniques while moving gradually to more complex models [6]. The detailed 
motivation and an account of the relations to the M-theory can be found there. 

One of the characteristic property of the supersymmetric quantum mechanics with 
the Yang-Mills potential is appearance of the continuous spectrum of non-localized states 
together with discrete, localized bound states [7]. The supersymmetric vacuum is believed 
to be in the continuous sector, with discrete spectrum beginning at some nonzero energy 
in general. Interestingly, the D = 10 (and not less than 10) quantum mechanics has also 
a threshold bound state at zero energy, which agrees with the M-theory correspondence. 
Apart of the D = 2 [2], [8], these systems are not soluble and the overall picture just 
outlined has accumulated over the years of intense studies of particular issues [9] -[21]. 

In Ref. [6] we proposed to use the standard Hamiltonian formulation of quantum me- 
chanics. To this end we have constructed explicitly the (finite) basis of gauge invariant 
states and calculated algebraically matrix representations of a Hamiltonian and other rel- 
evant observables (e.g., supersymmetry generators). This done, we proceeded to compute 
numerically the complete spectrum, the energy eigenstates, identified their supersymmet- 
ric partners, computed Witten index, etc. The method has an intrinsic cutoff - the total 
number of allowed bosonic quanta rismax ■ Since our basis is the eigenbasis of the occupa- 
tion number operators, the cutoff is easy to implement. It is also gauge and rotationally 
invariant, hence it preserves these important symmetries. Since the size of the basis, i.e., 
the dimension of the cut Hilbert space, grows rapidly with Hb, convergence with the cut- 
off is the crucial question for this approach. It turns out that in all cases studied there 
(i.e., the Wess-Zumino quantum mechanics, and D — 2 and D — A, SU(2), SYMQM) 
many nontrivial results were reliably obtained before the number of basis vectors grew 
out of control [22]. The approach applies as well to bosons and fermions being entirely 
insensitive to the notorious sign problem which plagues any Monte Carlo attempts to 
attack these systems. Later on the new, recursive method of calculating matrix elements 
significantly improved the precision of the solution of the D — 2 SYMQM and eventually 
inspired the exact, analytic calculation of the restricted Witten index for this model [23]. 
To make further progress one has to deal with the rapidly growing number of states. Of 
course this problem is most severe in the D = 10 model where some preliminary results 
for pure Yang-Mills system were nevertheless already obtained confirming for example the 
SO (9) invariance [24]. 

It this paper we have abandoned the brute force construction and diagonalization of 
the Hamiltonian in the whole (cut) Hilbert space. Instead, we have exploited fully the 
rotational invariance solving the problem separately for each angular momentum. Second, 
the recursive construction of matrix elements of Ref. [23] was generalized and adapted to 
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the fixed angular momentum channels (Section III). The two tricks coupled together led 
to the quantitative improvement of the precision and allowed to uncover a rich structure 
of the system to a much more complete level (Section IV). 

Finally, for the scalar (j = 0) sector, one can push the cutoff even higher performing 
complete analytic separation of variables in this case [25, 26]. Using the method adapted 
by van Baal for the noncompact system considered here, one can reach yet higher cutoffs 
in the two {np = 0, 2) channels. Results of this approach will be briefly discussed in the 
next Section. 

Recently, a new possibility to optimize numerical solutions for the lowest state of the 
system has been investigated [27]. 

Effective Lagrangians for various dimcnsionally reduced supersymmetric Yang-Mills 
theories, including SYMQM, have been very recently derived in Ref. [28]. 

Supersymmetric Yang-Mills theories in extended space have been studied for some 
time with the aid of the Hamiltonian approach on the hght cone [29]. 



2 The system and early results 
2.1 Definitions 

The reduced quantum-mechanical Yang-Mills system is described by nine canonically 
conjugate pairs of bosonic coordinates and momenta x^(t), Pait), i = 1, 2, 3, a = 1, 2, 3 and 
six independent fermionic coordinates composing a Majorana spinor 'i/'"(t), a = 1,...,4, 
a = 1, 2, 3 satisfying canonical anticommutation relations. In D = 4, it is equally possible 
to impose the Weyl condition and work with Weyl spinors. The Hamiltonian reads [30] 

H — Hk + Hp -\- Hp, 

. . . . 

Hp — —eabc^ade^b-^c-^d-^ei (-'-) 
Hp = '-^eabc^a^''^bX^c; 

in D = 4, r'^ are the standard Dirac a'^ matrices. In all explicit calculations we use the 

Majorana representation of Ref. [31]. 

Even though three-dimensional space was reduced to a single point, the system still has 
an internal Spin(3) rotational symmetry, inherited from the original theory, and generated 
by the angular momentum 

f = e^^' + ^VjS^'^a) , (2) 

with 

E^''^ = -^[r^r'^]. (3) 

Furthermore, the model posesses the residual of the local gauge transformation generated 

by 

Ga = eabc (x^Pc - ^^l^c^ , (4) 
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and it is invariant under the supersymmetry transformations with the generators 

The bosonic potential Hp in Eq. (1), when written in a vector notation in the color space, 
has a form 

V = X x^f, (6) 

which exhibits the famous flat directions responsible for a rich structure of the spectrum. 

2.2 Creation and annihilation operators 

The Hamiltonian (1) is polynomial in momenta and coordinates. Therefore it is convenient 
to employ the eigenbasis of the occupation number operators associated with all degrees 
of freedom. To this end we rewrite bosonic and fermionic variables in terms of creation 
and annihilation operators of simple, normalized harmonic oscillators 

K,«r]-e {firj} = K:, P,^=i,2, (7) 

obeying the canonical (anti)commutation relations 

[4,Pc] = ^^^'^(.c, {^t. V^f } = (8) 
As usual bosonic variables are given by 

4 = -^(aUa?), K = ^(a^,-a?). (9) 

For fermionic variables we use the following representation for a quantum Hermitian 
Major ana spinor 

~fa ~^ '''fa ~^ "^fa ~ fa 
\-^f'a-f! + f'J + ^f'JJ 

which is easily shown to satisfy Eq. (8), with the help of Eq. (7). Other choices of fermionic 
creation and annihilation operators are possible [13, 30, 15]. 

2.3 The cutoff 

For completeness, we shortly review the practical construction of the cut Fock space used 
in Ref. [6]. The entire Hilbert space is generated by all independent polynomials of the 
elementary creation operators a^^ and f^^ acting on the empty state, i.e., the state with 
zero occupation number for all of the above-defined oscillators. In practical applications 
we shall work in the finite-dimensional Hilbert space of states containing a total of at 
most ub max bosouic quanta, i.e., 

nB<nBmax, rig = ^ a'Jo'f,. (11) 

i,b 

There is no need to cut the fermionic number, which is limited to 6 by construction. 



/ 1 + ^ 



(10) 
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The physical Hilbcrt space is restricted to gauge-invariant states only. It can be 
conveniently generated by all independent polynomials of gauge-invariant creators - 
bilinear or trilinear combinations of a^^s and /^'s (the explicit form will be given later). 

Finally, since elementary creation and annihilation operators have a straightforward 
action in the occupation-number basis, one can readily calculate all matrix elements of 
the Hamiltonian and other obscrvables. 

All these steps can be implemented automatically in a computer algebra system. The 
matrix elements of any operator are calculated by writing the operator in terms of creation 
and annihilation operators. Finally, the complete spectrum and eigenstates of the cut 
Hamiltonian (1) are obtained by numerical diagonalization. 

This approach has proved reasonably successful. But of curse there is a limit to it. It 
is possible to improve the results considerably by exploiting fully the symmetries of the 
cut system, and by foregoing an explicit construction of states using only matrix elements 
(cf. Ref. [23]). 

2.4 The symmetries 

Some of the symmetries were already exploited earlier to reduce the size of the bases. We 
now discuss shortly their significance. 
The fermion number np is conserved: 

[nF,H] = 0, nF = r^ra. (12) 

This is best seen in the Weyl representation of Dirac matrices, where the Majorana spinor 
10 assumes the simple form [32] 

i^l-{fl-flhJ2). (13) 

Since the Dirac a matrices arc block-diagonal in this representation, the fermionic Hamil- 
tonian Hp contains only bilincars of the type /^/. Therefore it cannot change rip- Because 
the Pauli principle allows only six colored Majorana fermions in this system, the whole 
Hilbert space splits into seven sectors, nj? = 0, 1, . . . , 6. The cutoff on the bosonic quanta 
preserves np^ and consequently the diagonalization described above can be carried out 
independently in each fermionic sector for finite n^max- 

The system is invariant under the particle-hole symmetry 

riF <-> 6 — np, (14) 

therefore it suffices to find the spectrum only in the first four sectors, np = 0, 1, 2, 3, with 
the Up — ?> sector being selfdual with respect to Eq. (14). 

The local gauge invariance of the full (non-reduced) theory turns into a global con- 
straint of the reduced quantum mechanics. Namely, the physical Hilbert space consists of 
the gauge-invariant states, which in this case are invariant under the global SU(2) rota- 
tions in the color space. This is taken care of by using the gauge invariant combinations 
of the creation operators. This symmetry is preserved by the gauge invariant cutoff (11), 
and was already maximally exploited by working exclusively in the color-singlet sector. 

On the other hand, rotational invariance had not been fully used until now. Again, the 
cutoff is rotationally invariant and, accordingly, only exactly degenerate SO (3) multiplets 
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Table 1: Sizes of bases generated in each fermionic sector, np. Ng is the number of basis 
vectors with given number of bosonic quanta, ns, while E gives the cumulative size up 
to ub- The last column gives the difference between the total number of the bosonic and 
fermionic states in all seven sectors. 

with well-defined angular momentum were observed in the spectrum. However no attempt 
was made to generate separate bases in each angular momentum channel. This is the 
main source of improvement in the present work and will be discussed in detail in the 
next Section. 

The system is also invariant under parity. In the F — sector, it is equivalent to 
bosonic parity, (—1)"-^, and states can be classified according to the parity of ns- 

Finally, supersymmetry is broken by limiting ns, since the generators (5) do change 
the number of bosonic quanta. It is therefore interesting to look for the restoration of 
supersymmetry with the increase of the cutoff. Indeed this was qualitatively observed 
earlier. Present improvements reveal the supersymmetric spectrum with much better 
precision. 

2.5 Early results 

The pf structure of the fermionic Hamiltonian has an instructive consequence. The 
interaction term Hp vanishes in purely bosonic sector np = 0, which means that the 
effective Hamiltonian in this sector is just the pure Yang-Mills, zero-volume Hamiltonian 
which provides the starting point of the small volume expansion [33]. Indeed the lowest 
eigenenergies found in this sector agree very well with well established results of Ref . [34] . 
Later on, this test was extended to higher states crosschecking with recent results by van 
Baal to 15-digit precision [35, 36]. 

Sizes of bases which were reached in Refs. [6, 22] are quoted in Table 1. They contain 
all angular momenta up to jmax < ns + ^np, also shown in the Table. Due to the particle- 
hole symmetry the structure in the np = 4, 5, 6 sectors is identical with that in np = 2, 1, 
respectively. It will be interesting to compare Table 1 with our new results displayed in 
Table 2. 

In Fig. 1 we display the lowest eigenenergies as a function of the cutoff in all fermionic 
sectors. Clearly the cutoff dependence is different in the np = 0, 1 than in ni;' = 2 and 
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Figure 1: Spectrum of D = 4 SYMQM obtained in Ref. [22]. 



3 sectors. Based on the experience with simpler models, where the correlation between 
the nature of the spectrum and the rate of convergence with Ncut was established, it was 
claimed that the spectrum in the = 0, 1, 5, 6 sectors is discrete, while it is continuous in 
the "fermion rich" sectors with n^? = 2, 3 and 4. Recent analytic solutions of a sample of 
quantum mechanical problems in a cut Hilbert space have proven that indeed continuous 
spectra are characterized by the slow, power-like dependence on the cutoff [37]. All 
these early results provided an evidence that sizes of the bases displayed in Table 1 were 
sufficient to calculate lowest locahzed states with a reasonable precision. 

By computing directly supersymmetric images of lowest eigenstates it was found that 
SUSY in the cutoff system was broken on the level of 10 - 20 % . This was also confirmed 
by the Witten index calculation. 



2.6 Separation of variables 

The above conclusions, about the signature and coexistence of the discrete and continuous 
spectra, have been dramatically confirmed recently by van Baal [36]. Decomposing the 
solutions of the nine- dimensional Schrodinger equation, in the = 0, j = and np = 
2, J = channels, into covariant tensors, the problem was reduced to a numerically 
affordable set of coupled ordinary differential equations. As a consequence, van Baal 
was able to push a cutoff up to n^max = 39 in these two channels, as shown in Fig. 
2. The discrete, localized states with n^r = can be clearly seen with a very high 
precision. Moreover, the intricate nature of the np = '2 sector is also evident. As expected, 
the localized states have quickly convergent eigenenergies while the continuous spectrum 
manifests itself as a family of levels which slowly fall with the cutoff. We postpone the 
detailed discussion of this beautiful result until the global picture of the solutions in all 
channels becomes clear. 

Let us move now to the main subject of this paper which extends the results just 
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Figure 2: High cutoff results from van Baal approach for two available channels: np — 
2,j = (upper) and np = 0, j = (lower). The dashed line for = is the only 
odd-parity level in this energy range. 

presented. The new method allows to reach cutoffs in the range 18 < n^max < 23 in 
all fermionic sectors and for all angular momenta, providing at the same time detailed 
information on the supersymmetric interrelations between eigenstates. 
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3 The new approach 



We first present the basic features of the new algorithm, which allows us to push the 
computation much further. 

Rotational symmetry is exploited fully: all the objects in the computation, beside 
being gauge singlets, belong to irreducible representation of the rotation group Spin(3); 
this allows heavy use of the traditional machinery of Clebsch-Gordan coefiicients and 3j 
and 6j symbols. (In the following, several formulae will be used; they are reported in 
Appendix A). In addition, parity symmetry is used whenever possible. 

Vectors are never constructed explicitly; we build instead a recursive chain of identities 
between matrix elements of operators; this follows closely our algorithm for the D = 2 
case [23]. 



3.1 Gauge-invariant operators with definite angular momentum 
for the bosonic sector 

To avoid possible confusion, let us rename the bosonic creation and annihilation operators 
defined in Eq. (7) and dl respectively. In order to create states of a definite total angular 
momentum — j{j + 1) and = m, take the combinations 

1 1 

so that Ofe^^lO) is a state of angular momentum (l,m); now define — (— l)-'^+"*a^"*: the 
new creation and annihilation operators satisfy the canonical commutation rules 

[ar, arl = ^"''"''Sbc, [ar, «rl = V3c'^z,,s,,, (le) 

aj^^ and a™ transform as spin-1 triplets under rotations; they have odd parity (— 1)"^. 
(Here and in the following, ^ denotes the usual Hermitian conjugation applied to a single 
Jz component of an operator; e.g., al^ is the Hermitian conjugate of a^.) 

Prom it is possible to build the bihnear gauge- invariant operators a^^a^^, which 
are then decomposed in components of given angular momentum Aj^^; let us introduce 
the notation 

{■^hi ^h) i,m — ^ ^mim2m-^ji,miSj2^rn2J (17) 
mi,m2 

where R and S are arbitrary operators with definite rotational properties; Eq. (17) implies 
We can now define 

Am = hm = (-1)^+"^ A,',-^ = («,, S,) (19) 

where Aj „j is the Hermitian conjugate of Aj^m- Since A is a symmetric combination of 
a's, it has no j=l components, but only 1 j=0 component and 5 j=2 components; A\^^ 
and ^2,m transform as spin-2 quintets under rotations. 
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In order to express the commutation rules between A and A'^ , it is necessary to intro- 
duce the gauge-invariant "mixed" operators 

Bj,m = -m^, aJ),-,„. + (al a,),g = (-l)^^-^^; (20) 

in addition to 1 j=0 component and 5 j=2 components, B has also 3 j=l components. 
We can now write 

j,m 

[Bji,m2T ■Bj2,m2\ — ^^^jij2j^mim2m-Bj,m- (23) 



It would be pointless to write the detailed form of the coefficients c, c', and c"; their 
computation will be discussed in Appendix C. 

We also introduce the trilinear gauge-invariant creation and annihilation operators 

= Sbcdal^al^al^ = i£bcd{{al, al)^, al)^^, A^A (24) 

(the notation ((-Rj^, Sj.^)j.^,Tj^)j.rn follows from applying Eq. (17) twice), which have only 
the scalar (i.e., spin-0) component, and the "mixed" trilinear operators 

^J,m = £6cd((oJ, 4)^, a<i)^. ^, Bj^rn = Sbcdiih, ac)i, 4)j,m' ^hrn = (-l)^^'"^^-,™. 

(25) 

The above-defined operators form a complete set of gauge-invariant bosonic operators, 
in the sense that any gauge-invariant bosonic operator can be written as a polynomial in 
these operators. In particular, we can write Hk and Hp in terms of A} , A, and B as 

Hk = \V2>{Alo + Ao,o - 25o,o); (26) 
j^Hp = \{A\ At)p + i(i, A)p + (At, B)p + (S, A)p + (S, S)p + |(At, A)p 

+ V3(4,o + ^,o + 5o,o), (27) 
{R, S)p = 2i?o,0'S'o,o — \/5(-R2, 'S'2)o,o- 



3.2 Fermionic operators with definite angular momentum 

To identify fermionic creation operators with definite angular momentum, recall the origin 
of the parametrization (10). It represents a Majorana fermion in Majorana representation 
of Dirac matrices and was obtained by a unitary transformation of a Majorana fermion 
in the Weyl representation (13) [6]. Therefore /^t creates in fact a fermion in the Weyl 
representation and as such carries definite angular momentum. This follows from the 
explicit form of the spin operator 5*^ defined in Eq. (2): 

= \€^''i^a = \ [fVfl - fl'ft) , (28) 
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which can be obtained in either Weyl or Majorana representations of Dirac matrices. 
Therefore, the fermionic creation and annihilation operators defined in Eq. (7), are 
already the desired operators, and we set 



1 ^ 1 ^ 

?1 ^"2 _ ^2. 



f2 fi f 2 f^. 

Jb — Jbi Jb — Jb-i 

and /™ are j = \, Jz = Tn operators; the anticommutation relations are 

f fmi fm2t\ _ X X, r fmi fm2U__/n/^hh^ r 

\Jb tJc J — "mim2"6c) \J b tJc J— ^ ^mim20 ^bc- 

3.3 Gauge-invariant operators involving fermions 

Let us complete the set of gauge-invariant operators, bihnear C^^'"^ or trilinear O^^^'^ , 
with definite J^, J^, n^, and ub- In the bilinear case, we complement the bosonic opera- 
tors 0]^l = A\^, a, and = Bi^ with 

o^'-^ = e^ = 

(G'im vanishes identically); note that F, G, ©, and $ give zero when applied to a bosonic 
state. 

In the trilinear case, we complement the bosonic creation and annihilation operators 
Oo'o = and A with 
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the antisymmetrized product of two a"^'s only produces j = I, and likewise for /''"'s; the 
factor i in F^, (5^, and P is required to have real matrix elements of H between ans even. 
Up — 1 and a riB odd, n^? = 1 state. We also define the "mixed" operators — B\ B, 
and 

^j,m, ~ ^i,m ~ ~'^^bcd{{fbjfc)lifb)j,m^ ^j,m = ^j,-m) J = 2' 2' 

0?^' = = -^Bbcdiifl fc)f, = / = 0' 1' J = 0' 1' 2; 

Ojlm — "],m = ~i^bcd{{fbJ fl)lj ^d)j^rn^ ^j,m — { — ^y~^"^^j,-m, j — 0, 1, 2; 

O]^' = = -i^bcd((a6, 4),-,, /]). ,„, Irj'^ = (-l)^+-^r/_^, f ^ 0, 1, 2, j = i, I, f ; 

Oj,m ~ ®],m ~ "i^bcdiiflbiOjcjl-, fd) j,mi ^3,m — ("I)"'"'' ^j,-mi j = 2' 2' 

note that Q^, l», S, ^, and © give zero when apphed to a bosonic state. 
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We can establish (anti) commutation relations between pairs of gauge-invariant oper- 
ators, similar to Eqs. (21) and (22); it would be pointless to present here their explicit 
form; their computation will be discussed in Appendix C. 

The above-defined operators form a complete set of gauge-invariant operators, in the 
sense that any gauge-invariant operator can be written as a polynomial in these operators. 
In particular, we show by explicit computation that 

HF^-gVSiHl + ^io)- (29) 

Q must be decomposed in components with definite and np', since in our Majorana 
representation P is real and T,^'' is purely imaginary, Qa is Hermitian. Denoting the 
Up — i and Up — —i doublets by Ql^ and Qm respectively, we have 

Q\ = le''{Ql-^Q2 + Q3 + ^Q4), 

2 

Qm = (-l)^+'"g_„, (30) 
where 9 is an arbitrary phase; the anticommutation relations are 

{QLQl} = 0, {Qm,Qm'} = 0, 

where is defined by the analogous of Eq. (15); f „j gives zero when applied to a gauge- 
invariant state; the only nontrivial anticommutator can be rewritten as 

H^lJ2{Qm:Qlh (31) 

m 

we choose 9 — — ^tt; an explicit computation gives 

q1^ = 7|(fJ -el -IgFt +i^et +gW). (32) 
Note that, with the present conventions, all matrix elements of interest are real. 

3.4 Construction and orthonormalization of states with definite 
angular momentum 

All states are classified into even and odd states, according to the parity of p = ni? -|- 
ub (mod 2). (This label coincides with parity only for np — ^ states.) 

It is useful to set up a common naming scheme for all our creation operators: X{y^-p)'^ 
is the creation operator with np = v and = 2 + p — z/; i.e., X(0, 0) = A, X(l, 0) = F, 
X(2,0) = G, X(0,1) = A, X{\\) = F, X(2,l) = G, and A:(3, 1) = /; A:(3,0) is 
identically zero. 

We build our states recursively, applying Xiy^]))^ to a state of an orthonormal basis 
with definite and J^, and taking linear combinations to produce again an orthonormal 
basis. 
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It is important to note that, given the e contraction rule 



^iii^h ^hhh ~ ll'^iaj/all) (33) 

the product of two trihnear operators can always be decomposed into a sum of products 
of three bilinear operators; therefore, even states can be built by applying any number of 
even (p = 0) creation operators to the vacuum; correspondingly, odd states can be built 
by applying one odd (p = 1) and any number of even creation operators to the vacuum. 

is never needed in combination with fermionic operators, since X(z/, 0)^74^ for > 
can be written as a linear combination of terms of the form IjMt. 

Creation operators (which (anti) commute between themselves) can be ordered to get 
every fermionic operator to the left of every bosonic operator and every trilinear oper- 
ator to the left of every bilinear operator. Therefore, using the notation for 
our states, we can build all bosonic states from even bosonic states as |0,2n+p) = 
X(0,p)t|0, 2r2— 2), and all fermionic states of parity p from even states of lower np as 
= X{u,p)'^\np—u,nB—'2,—p+i'), 1 < < rip. 

In order to create a fermionic state with ub > np -\- 1, sA, least one must be used; 
therefore, such states can also be built as \nF,nB) — A'^\np,nB—'^)\ this second recipe 
turns out to be much more efficient, both in generating and orthonormalizing the states 
and in computing matrix elements of operators between them. 

A basis for the sector with given np and is contained in the set 



\j,m,np,nB;i^,p,ji,j2,i)^ J2 C^m^mimX{^^P)h,mi\h,'rn2,np-u,nB-2-p+u;i), (34) 



rr, 

m\,m2 

where 



i/ = 0, p = np -\- ub (mod 2), np — 0; 
1 < i/ < np, p = np + nB (mod2), np > 0, nB < np + 1; (35) 

u — 0, p — 0, np > 0, nB > np + 1. 

The scalar product of two such states can be written as 

if, m', n'p, u'b] ly', p', j[j2, m, np, n^; i/, p, ji, ja, i) = 

Oj'jOm'mOn'pnFOn'j^nB '^u' ,p' ,j[,j^,i'■,u,p,jl,j2,^■ 
By Gram-Schmidt orthonormalization we obtain the orthonormal basis 

\j,m,np,nB;i) = ^ Ri:^^;j,%,j\j,m,np,nB;iy,p,3i,j2,j), (36) 

Epj,nF,nB qj,nF,nB r>j,nF,nB r 

i^' ,P' ,j[ ,32 ' ^'P'jl '3 

The states of the set (34) may not be linearly independent; this is however not a serious 
problem: Gram- Schmidt orthonormalization will select an orthonormal basis and give a 
non-square matrix R. Eq. (36) implies 

(j, np, nB-, i'\\X{u,p)]^ || js, np-u, nB-2-p+u; i) 

^ ^ 2^ ^i'y,p\j[j!,j''^u',p\j[,j'^,i'-,u,p,n,32,v K'^') 

t^' ,P' ,3[,3'2'i' 
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where denotes a reduced matrix element, cf. Appendix B. 

To compute the scalar product of two such states, define 

where the sign is + (anticommutator) when both u and i/' are odd, — (commutator) 
otherwise. 

Using completeness and applying well-known identities similar to Eqs. (48) and (55), 
we obtain 

(/, m', np, ub] i^',p',j'i,j2, m, np, ub; '^,P,ji,j2, i) 

j3,i3 U3J1J2J 

X {j2, rip, n'B] i'\\X{v,p)]jj^, n'p, n'^; is) {j^, n'p, n'^; i3\\X(u' ,p')j^ \\j2, np, ub; i) 
+ 5,,5^,^Y.i-ir''^^^^'''i^l I I } 02,ri^,r^^;^1l<'''"^^'^''''■^^llJ2,n^,nB;^), (38) 

where rip, ub, n'p, Ub, 'n'p, and n'^ are fixed by the selection rules. The r.h.s. involves 
matrix elements of operators between states with lower np or ub- 

In the case p — p' — 1, the (anti)commutators would involve all trilinear operators; to 
avoid this, applying Eq. (33) we rewrite X^u', 1) X(z/, 1)''' as a sum of products of three 
bilinear operators, decomposed in components of definite angular momentum; they are 
dealt with exactly like the commutator term in the above equation, with the same factors 
and 6j symbols. The explicit computation of the decomposition will be discussed in 
Appendix C. 

3.5 Recursive computation of matrix elements of operators 

Our task is to compute a matrix element of the form 

{f, m', n'p, h'b; i'\Oj"^rn"\j, m, np, nB] i), 

where O is an operator with a definite number of fermionic and bosonic quanta n"p and 
n'^; let us take n'p = np + n'p, n'^ = ub + n'B (otherwise, the matrix element is zero). 
Apply Eqs. (34) and (36) to the ket: 



(/, m', n'p, n's; i'\Oj»^rn"\j, ni, np, nB] i) 

Efjij2j r>j,nF,nB 

X (/, m', n'p, n'B] i'\Ojn^^„X{v,p)] jja, ^2, np-v, nB-2-p+u; j) . 



L nL- i'\n.„ „ Y(ij Ti, . 

Using the (anti)commutator 



33,1^3 
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and completeness, we obtain 



(/, n'p, n's; i'\\Of'\\j, uf^ub] i) 

!',pjij2j;j3,»3 ^•''^ •' •' > 

X (j', ^b; p)]^ II j3, n^-z/, n'5-2-p+i/; ig) 

X (is, rip-i^, ns-2-p+i/; isljC^v/jj, j2, np-v, nB-^-p+u; j) 

X (j', n^, n'^; ^'ll^jg^^''''^'^'^ 11-^2, n^-i^, nB-2-p+z/; j). 

In the case of p = 1 and trilinear (9, we can again resort to the use of Eq. (33) to rewrite 
O X{v, 1)^ as a sum of products of three bihnear operators, decomposed in components 
of definite angular momentum. Every matrix element is computed in terms of matrix 
elements for smaller np and/or n^; the recursion is closed when a matrix element is 
obviously zero, or when Eq. (37) can be applied; the only nontrivial case is 

(0, 0, 0; 1||5^||0, 0, 0; 1) = -\V^6jQ. 

Finally, to compute H apply Eqs. (26), (27), and (29); to compute Q, apply Eq. (32). 
The implementation of the algorithm will be described in Appendix C. 



4 Results 

4.1 Hilbert space: sectors, channels and diamonds 

The approach described in Sect. 3 allows to deal with a considerably larger Hilbert space 
than the direct method, cf. Tables 1 and 2. With the recursive algorithm implemented in 
Mathematica we were able to compute all matrix elements of H and Q on a single PC in 
a time ranging from 2 minutes for np — ^ alone to 140 hours for the whole computation. 
The whole Hilbert space was effectively spht into seven sectors of fixed fermion number, 
< < 6, which in turn decouple into channels of fixed angular momentum j. In 
Appendix D we quote sizes of bases in all {np,j) channels for all available values of Ub- 
It is useful to represent this decomposition on a {np,j) plane where circles corre- 
sponding to the individual channels form a regular mesh ^. Fig. 3 shows such a map 
together with the number S0(3) multiplets in each channel (for a particular value of a 
cutoff Ub max )• The distribution of states among channels is such that each SO (3) mul- 
tiplet belongs to one and only one diamond adjacent to the vertex. This "population" 
of all vertices determines the multiplicities of individual diamonds, i.e., the number of 
supermultiplets reproduced at given ub max ■ The latter are given in italic in the figure. 
More precisely, if Rj denotes the multiplicity of a diamond I, and di is a number of S0(3) 
multiplets in a channel i — {np,j) then 

d, = T^^iRi, (39) 
^With two symmetric dislocations: there are no states with j = 1 and np = 0,6. 
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Table 2: Number of S0(3) multiplets (left) and size of the basis (right) of the physical 
Hilbert space, at given ns for each np sector. Right columns correspond to the left ones 
of Table 1. The last column gives respective sizes summed over all seven fermionic sectors. 

where S/|j means summation over I's adjacent to i. These diamonds are nothing but 
super symmetric multiplets which will be discussed below in detail. At present stage 
they play only a kinematical role — they provide an alternative way of classifying all 
basis states. Our cutoff riBmax violates supersymmetry, however it violates it rather 
gently. Namely, for every odd nsmax the total number of states is such that they fill 
the integer number of supermultiplets. Taking into account the SO (3) degeneracy, we 
see immediately that the numbers of fermionic and bosonic states in a diamond match. 
These last two properties account for the exact balance between fermionic and bosonic 
states found earlier, cf. the last column of Table 1. For even riBmax j Eq- (39) also holds, 
but exactly two diamonds at highest j have d = ~1. However, this does not spoil the 
exact balance, and it is irrelevant in our perspective of increasing nBmax at fixed j. 

We shall discuss now some detailed features of the model which follow from supersym- 
metry and rotational symmetry. 
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1 2 3 4 5 6 np 

Figure 3: Populations of channels (bold) versus equivalent number of supermultiplets 
(italics) at ub^,, = 3. 

4.2 The algebra of SUSY generators: supermultiplets 

It is convenient to work with the Weyl generators defined in Eq. (30), which carry def- 
inite fermionic number np and angular momentum (m = ±|), and satisfy the standard 
anticommutation relations 

{Qm,Qi}^^HSrn^, {Ql,Qi}^{Qm,Qn}^0 (40) 

in the gauge-invariant sector. 

A supermultiplet can be constructed by considering Ql^ as creation operators and Qm 
as annihilation operators acting on the "vacuum" - the lowest-n^r state of the supermul- 
tiplet. Starting from a single eigenstate of H with nonzero energy and definite ni? < 3, 
j, and m, Eq. (40) implies that exactly three more states are produced: two by acting 
with Q\ and i , with quantum numbers np + 1 and m ± |, and one by acting with 

2 ^2 

Q\ 1 1 Q\ ) with quantum numbers np + 1, j, and and m (since Q\ i i Q\ 

2 ~2 ~2 2 2 ~2 ~2 2 

vanishes). It is easy to show that applying or Ql^ to any of these states either gives 
zero or another of these states. Starting now from a full rotational multiplet, \nf,j, m) 
is decomposed into two multiplets with np + I and j ± |, (j — | is absent if j = 0); addi- 
tionally, we have one np,j and one np + 2,j multiplet. This structure, which is nothing 
but a diamond of Sect. 4.1, is shown in Fig. 4. 
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Figure 4: The diamond structure of a supermultiplet, with supersymmetry fractions q. 



4.3 Supersymmetry fractions 

Applying Eqs. (57) and (58) to Eq. (40), we obtain 

^TT^i^E [{r,iV\\f;i') {r,imj';i') + -(41) 

When \j;i) is an eigenstate of H with eigenvalue Ej^i and i" — i, Eq. (41) reduces to 



E 



i'- .-'lint I 



= 4(2j + 1)£;,- 



(42) 



It is therefore natural to define the "supersymmetry fractions" 



{f;^'\\Q^\J;^) 



which, in the limit of exact supersymmetry, satisfy the sum rule 

Efe(/'^1j'0 + ?a^l/,O] = 2j + i. 



(43) 



Supersymmetry fractions allow an easy classification of states into supermultiplcts: 
q is nonzero only between states belonging to the same supermultiplet. For discrete 
states, the fractions can be computed explicitly, using Eq. (42) and remembering that 
OVII(O) vanishes on gauge-invariant states. When no mixing occurs, the re- 

sulting fractions are 



q{nF+l,j+^\nF,j) = q{nF+2, j\nF+l, = j + 1, 
q{nF+lJ-^\nFj) = q{nF+2J\nF+l, j-^) = J, 



(44) 



(they are also shown in Fig. 4); they saturate Eq. (43). 
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Figure 5: Cutoff dependence of the spectrum in a sample of channels. 

4.4 The spectrum - cutoff dependence and general properties 

Monitoring the cutoff dependence is crucial for at least two reasons. First, it provides a 
model-independent information about the errors induced by limiting the Hilbert space. 
Second, it allows to distinguish between localized and non-localized states. The latter 
feature is particularly useful in studying supersymmetric gauge systems where continuous 
and discrete spectra are known to coexist. It was shown in Ref. [37] that eigenenergies 
of non-localized states drop slowly to zero with the cutoff, while the discrete spectrum is 
characterized by rapid convergence to the finite, "infinite volume" eigenvalues. 

In the (0, 0) and (2, 0) channels our results are identical with those of van Baal, Fig. 2, 
hence we concentrate on other channels, plotted in Fig. 5. For nj? = 1, j = | the spectrum 
is similar to that in the (0, 0) channel. The levels are quickly converging and the available 
cutoff is sufficient to guarantee small errors. Similar situation occurs for higher angular 
momenta with np = 1, e.g., for J = |. Some degeneracies are observed, e.g., second and 
third level of the channel {np = l,j = |). They are not caused by supersymmetry, which 
connects states from different channels, as was discussed in Sect. 4.2. 

On the other hand in the np — 3 sector we clearly observe both convergent, localized 
states and slowly falling ones from the continuum. The j — \ channel plotted in Fig. 5c is 
very similar in this respect to the np = 2, j = channel shown earlier (cf. Fig. 2). Similar 
behavior is seen for other angular momenta. Again, cutoffs reached with present method 
allow for quantitative studies of many features of the locahzed states. Scattering states 
show much more complexity, nevertheless some of their properties can be also inferred, 
see below. 

We therefore seem to confirm the general pattern suggested already by the low n^max 
calculations: in zero- and one-fermion sectors (and their particle-hole images) the spec- 
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Figure 6: Spectrum of the j = 0, ^ states obtained with highest cutoff available in each 
fermionic sector, cf. Table 2. 

trum is discrete, while in the "fermion rich" sectors with np = 2,3,4, both localized and 
non- localized states coexist. There is however additional refinement of this rule. 

Contrary to earlier expectations the spectrum is entirely discrete also in the np — 
2,j = 1 channel, Fig. 5d. In fact we observe that this happens in all np = 2 channels 
with odd angular momentum j. Therefore previous rule is modified to the following: 
scattering states exist in the = 3 sector for all angular momenta, while non- localized 
states with np — 2 occur only for even angular momentum. This will find yet simpler 
interpretation when we discuss in detail the supermultiplet structure of the spectrum. 

4.5 Discrete spectrum - identifying supermultiplets 

To begin with, let us collect the "spectroscopy" graph of the energy levels with lowest 
angular momentum in all fermionic sectors. Fig. 6. 

Clearly a number of states in adjacent channels have identical energies (within our 
cutoff errors) and are therefore good candidates for SUSY partners. Confronting this 
with the cutoff dependence. Figs. 2, 5, we see that identification of SUSY multiplets is 
simpler in the discrete part of the spectrum. Restoration of supersymmetry among the 
non-localized states is more complex and will be discussed later. Still, in order to achieve 
a complete classification (even of localized states), it is important to analyze together the 
highest cutoff results, the cutoff dependence, and supersymmetry fractions. This is done 
below. 

Recall from Sect. 4.2 that a supermultiplet of SYMQM is composed by the diamond of 
0(3) multiplets shown in Fig. 4: {np,j), the multiplet with the lowest np, {np+1, 
{np+1, j — ^) (only if j > 0), and {np+2,j). We will denote the full supermultiplet by 
the spectroscopic labels np{j) ^ for the ground state in the channel, np{j)' for the first 

^Note that the pairs of supermultiplets conjugated by particle-hole symmetry are 0{j) with 4(ji) and 
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Figure 7: Energy levels vs. ub max for the channel (0, 0); solid lines indicate even n^-parity, 
dashed lines odd n^-parity. 

excited state etc.; when many excited states are considered, we label them by their energies 
multiphed by 10^. In the case np — 0, when (— 1)"'^ is conserved for the (0, j) multiplet, 
we add an n^-parity label, i.e., we write 0{j'^). 

We begin the detailed presentation of our data by plotting the energy levels vs. n^max 
for several channels in Figs. 7-11. The channels with higher j follow a pattern quite 
similar to these, and therefore we will not present here the corresponding plots; for the 
remaining channels with j < 4, they can be found on the authors' web site [38]. For 
the lower levels of each channel, we quote the spectroscopic labels npiJ) identifying the 
super multiplet to which the level belongs, anticipating results from the following of the 
present Section. 

The most effective tool to classify states into supermultiplets is based on the analysis of 
supersymmetry fractions; the spectroscopic labels reported in the above-mentioned plots 
are obtained by the following method. 

Let us select two sectors with fixed {n'p = nF+l,j' = j±^) and {nF,j), and construct 
the matrix q{n'p, j' ,i'\nF, where i' and i run over the energy eigenvalues of the two 
sectors; the cutoff ub max is often the same in the two sector, but it may be different, in 
which case we will write the two cutoffs as n'^^g^^\nB max- Take one (n^, /, i') state and look 
at the corresponding row of the matrix as n^niax grows: if all elements go to zero, the state 
belongs to a {n'p,j') supermultiplet. In the same way, take one {np,j,i) state and look at 
the corresponding column of the matrix: if all elements go to zero, the state belongs to a 
(n^— 2,j) supermultiplet. If one element remains nonzero, the two corresponding states 
belong to the same supermultiplet: we look for the remaining supcrpartners coupling to 
these two states in the appropriate channels, forming the diamond of Fig. 4, with the given 
values of q. If two elements remain nonzero, we have a case of "accidental" degeneracy 
of two supermultiplets: the g's are the superposition of two patterns of Fig. 4, with 

with 3(j), while 2{j) is self-conjugated. 
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Figure 9: Energy levels vs. n^max for the channel (2, 0). 

coefficients cos^ 9 and sin^ 9, where 9 is the mixing angle between the energy eigenstates 
(which are not exactly degenerate at finite ub max) ^-nd the states belonging to a definite 
supermultiplet. If a number of elements remain nonzero (typically 5 to 10 for our values 
of nsmsix), the state belongs to the continuum. 

Let us look in details, e.g., at the transition g(l, ||0, 0"*"). The q matrix for our highest 
value of TiBmax IS showu in Table 3 for n^max = 18; selected coefficients are plotted vs. 
nBmax in Fig- 12. We identify the states in each channels by their energies at riBmax — 18 
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Figure 10: Energy levels vs. rig max for the channel (2, 1). 
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B max 



Figure 11: Energy levels vs. nsmax for the channel (3, |). 

(jT-Bmax = 19 for (0, multiplied by 10^: we use the notation {np,j, E), or just E when 
np and j are obvious. 

Proceeding by increasing energies, first we see a perfect match for the two ground 
states, i.e., ^(411714117) = 1, and they therefore belong to the O(O^) super multiplet. 
Next we have a case of mixing: since g(6388|6388) = 0.114 and g(6401|6388) = 0.885, 
(0, 0+, 6388) belongs to 0(0+)', while (1, |, 6388) and (1, i, 6401) are linear combination of 
states belonging to the "accidentally" degenerate 0(0+)' and (1, |) super multiplets, with 
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Table 3: Matrix of supersymmetry fractions 5(1,110,0+) at UBmax — 18. States are 
identified by their energies. All numbers are multipled by 10^. 
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Figure 12: Selected supersymmetry fractions g(l, |,i'|0,0+,i) vs. ub 
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Table 4: Matrix of supersymmetry fractions g(2, 0|1, |) at nemax = 18. 

a mixing angle ^ with cos^ ^ = 0.114. For higher levels, g's are not completely stable in 
^Smax, and we need to extrapolate them to nsmax oo. We clearly see that (0, 0+, 7997) 
and (1, |, 8063) belong to 0(0"*")". For higher states, the analysis requires more care. 

The remaining members of the 0(0"'") supermultiplets can be identified by looking at 
g(2, 0, i'|l, |, i) , which is shown in Table 4 for n^max = 17|18: (2,0,4121) belongs to 
0(0+); since g(6397|6404) ^ and g(6484|6388) ^ 0, (2,0,6397) and (2,0,6404) are 
linear combination of states belonging 0(0+)' and (1, |), with the same mixing angle 6 
as above. (2,0,8806), despite the high energy, is very easily attributed to 0(0~), with 
the help of Table 5. The levels (2, 0, i') related to the continuum spectrum have zero 
q{2,0,t'\l,l,t). 

We then look at g(2, 1, i'|l, |, i) and g(3, |, i'\2, 0, i) to identify the remaining members 
of the 1(|) supermultiplets. g(2, 1, i'| 1, |, i), presented in Table 6, presents the same 
pattern as g(2, 0, i'|l, |, i), except for the absence of continuum states, and we will not 
delve into the classification of states. 

g(3, |, ■i'|2, 0, i), shown in Table 7, presents a new, very interesting pattern: wc see 
states with a broad distribution of g's quite different from zero, even with states with 
very different energies; looking at the n^max dependence of the levels, cf. Figs. 9 and 11, 
we conclude that the patterns identifies continuum levels. On the other hand, q is zero 
between continuum and discrete states, or between discrete states of significantly different 
energies. We can easily identify members of supermultiplets with quantum numbers 
0(0), 1(|), 2(0), and 2(1), with the remaining states belonging to the continuum. The 
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Table 5: Matrix of supersymmetry fractions g(l, ||0, ) at nsmax = 18|19. 

"doubling" of (3, |) states belonging to 1(|) supermultiplets is due to the particle-hole 
symmetry, as will be explained below. 

It is also worth presenting g(3,|,i'|2,l,i), shown in Table 8; thanks to the absence of 
continuum states from (2, 1), it is very easy to identify states in (3, |) belonging to the 
supermultiplets 1(|) and 2(1). 

The analysis of q for higher j is repeated exactly in the same way. We will not present 
here the corresponding q matrices, which can be found in Ref. [38] for the remaining 
channels with j < 4. We only remark that all q{n'p, j) for the same values of n'p, np 
and different values of f,j are qualitatively very similar (in the case of = 2 (np — 2), 
only for j' (j) having the same parity). 

From all the above data, we can compile the spectroscopy of Tables 9 and 10. The 
table is limited to np < 2, since the other supermultiplets can be obtained by particle-hole 
reflection, and to j < 4, since nothing new happens for higher j. 

One feature should be stressed: for each supermultiplet, the particle-hole sym- 
metry implies the existence of a conjugate supermultiplet 3(j), and therefore of two (3, j) 
states of degenerate energy (in the riBmax limit); we observe mixing of each pair, 

with a mixing angle 6 = ir/A. 

Figure 13 shows a sample of the lowest supermultiplets for the first few angular mo- 
menta and all np- Degenerate supermultiplets at ~ 6.4 and 8.1 were shghtly spht for 
the sake of illustration. 

4.6 Continuous spectrum 

We already mentioned that a continuous spectrum is observed for all j's in the np — ^ 

channels, but only for even angular momenta in the np = 2,4 sectors. This pattern 
is consistent with supersymmetry, and simply means that continuous states exist only 
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Table 6: Matrix of supersymmetry fractions ^(2, 1|1, |) at ub max = 18. 

in supermultiplets 2(j) with even j. Note that the "opposite" behaviour (all np — 2 
channels and every second np = S channel) cannot be accommodated into a geometric 
structure of supermultiplets, cf. Fig. 3. 

It is also interesting to realize that, even though supersymmetry is broken by the 
cutoff, the above rule is not, i.e., we don't see any hint of continuum states in the channels 
{uf — 2,4, j — odd) for any finite cutoff. 

4.6.1 Scaling 

Non-locahzed states of the system describe D-particles [13] penetrating the fiat directions 
of the potential, as mentioned in the Introduction. In a cut system all energy levels of 
the continuum states fall to zero with increasing ub ma.x ■ If we label them by a "principal 
quantum number" m, the large-cutoff limit at fixed m is trivial. Such a phenomenon was 
also found in the free case when one regularizes the system by limiting the number of 
quanta [37]. In that case, it was also shown that the nontrivial and correct continuum 
limit is the scaling limit 

E{p) = lim m = ^V2N. (45) 
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Table 7: Matrix of supersymmetry fractions g(3, ||2, 0) at nsmax — 18. 
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Table 8: Matrix of supersymmetry fractions q{3, 1|2, 1) at n^max = 18. 
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Table 9: Spectroscopy of SYMQM. 

where p is the continuum momentum and the cutoff. These results were obtained 
analytically for a free particle in one dimension. They also apply to the D = 2 SU(2) 
SYMQM, since this is effectively a quantum mechanics of a free particle in three (color) 
dimensions, projected on the singlet and triplet channels of angular momentum [39] . The 
present, D = 4, case is more complicated. However, we expect that, whenever it is possible 
to define asymptotic states with given momentum, as is the case for the scattering process 
considered here, some version of Eq. (45) should hold. Scattering states in the present 
model correspond to particles propagating freely in the three dimensional (color) flat 
valleys of the potential V, cf. Eq. (6). Gauge invariance restricts color orbital angular 
momentum to few channels, so we are not that far from the D = 2 example. We have 
therefore taken Eq. (45) as a phenomenological rule and tested it with our data. 
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Table 10: Spectroscopy of SYMQM (continued). 

* The two states (2, 3, 8420) and (2, 3, 8527) belong to supermultiplets with different energies, but appear 
to be mixed at the available values of nsmax- 

The scaling limit (45) implies that at fixed m all energies of non-localized states behave 
as 0(l/iV). Figure 14 tests this prediction assuming that we identify the one dimensional 
cutoff N with riBmax ■ Indeed the energies of the first four levels seem to follow '^/nsmax 
behavior both vanp — 2 and np — sectors. We did not use higher levels since they are 
probably infiuenced by the the discrete spectrum. 

One can also contrast the m dependence with the one dimensional formula 

and with the D = 2 case. Table 11 compares ratios of our first four energy levels, for the 
largest value of the cutoff, with analogous ratios of the D = 2 system at the same value 
of A^, and with Eq. (46). 

The comparison is done in two channels: the (2, 0) channel, corresponding to the 
bosonic np — ^ sector of the D — 2 model, and the (3,1/2) channel, which is the 
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Figure 13: Spectrum of a sample of supermultiplets identified in Tables 9 and 10 
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Figure 14: Eigenenergies of the first three levels from the continuum, in the [np = 3,j = 
1/2) channel, multiplied by n^max , as the function of the cutoff. 

counterpart of the fermionic np = ^ sector. ^ To give an idea of the cutoff effects, we 
quote the D = 2 energies for = 18 = riBmax (third column) and for N = 150, which 
is easily available in this case and coincides with N = oo within the two digits accuracy 
reported (fourth column, lower half). 

■^See Ref. [23] for the details of the D — 2 system 
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Table 11: Ratios of the energies from the continuum. Comparisons between the D — 2 
and D — A systems. 

In the scalar case, high-cutoff results for D = 2 are identical with the exact ratios 
k'^/rn^. The onc-dimensional formula (46) does not apply to the fermionic sector. It is 
not surprising since this case corresponds to color angular momentum j = 1 and the three 
dimensional Schrodinger equation coincides with the one dimensional one only for j — 0. 

Finally, the comparison of the ratios for D = 4 with D — 2 is rather satisfactory. 
Numerical values of the energy ratios for the two systems are quite similar, over a range 
of an order of magnitude. All discrepancies are consistent with the cutoff effects. However 
one cannot exclude differences ~ 10% and consequently higher UBmax are required for 
more quantitative conclusions. 

4.6.2 Dispersion relation 

An interesting question appears whether the dispersion relation for the scattering states 
has the standard parabolic form, or whether it is modified by rather unusual behaviour 
of the potential. With the help of the scaling relation (45) we can now address this issue 
in both bosonic and fermionic sectors. For np — 2,j — the dispersion relation was first 
obtained by van Baal [35]. 

In Fig. 15 we have plotted the first three energy levels, as a function of rn/ y/n^^, 
for both bosonic {np — 2,j — 0) and fermionic {np — 3,j = 1/2) channels. Points 
from different m and n^max follow roughly a common curve which again confirms ap- 
proximately the scaling relation (45). Moreover, when the proper normalization of the 
momentum, required in (45), is taken into account, one obtains a reasonable agreement 
with the standard p^/2 kinetic energy of one degree of freedom (solid fines). 

Many effects prevent us from reaching better agreement at the moment. For example, 
the repulsion of the lowest discrete state at = 4.12 is clearly seen in the (2,0) channel, 
while it is not as efficient in (3,1/2), where the lowest state is higher [E — 5.19). For the 
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Figure 15: Dispersion relation for the lowest three scattering states in the (2, 0) and 
(3, 1/2) channels. 
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Figure 16: Regularized Witten index Iw{T) vs. T for even (left) and odd (right) nsmax < 
17 . 

present values of ub max , only the three lowest states can be used, hence one expects non- 
leading corrections in m. The identification of N with Ub max should be more carefully 
examined, etc. However, keeping in mind all these limitations, the overall picture seems 
reasonably satisfactory and we are looking forward for better data to make more extensive 
study of these points. 

4.7 Witten index 

With complete diagonalization of the Hamiltonian achieved in all sectors we can now 
calculate the regularized Witten index directly from the definition 

/^(T) = ^(-i)«Hi)e-Ti?W_ 

i 

The results, shown in Fig. 16, nicely confirm and strengthen early expectations based on 
much smaller n^max ■ 
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Figure 17: Infinite n^niax limit of the index deduced from the Fade approximants. 

As aheady mentioned, the number of bosonic and fermionic states is the same for any 
value of the cutoff, be it even or odd. Therefore the index vanishes at T = with this 
regularization. The sharp structure around T = clearly moves toward the origin with 
increasing cutoff indicating singularity at infinite ub max- Such a discontinuity is expected 
on general grounds and finds a reasonable support here. 

Of course at high T our "cut" index is bound to vanish exponentially. However there 
exist a range of intermediate times where definite fiattening occurs. This signals effective 
cancellations among supersymmetric partners hence a gradual, global restoration of SUSY. 
Moreover the plateau seems to converge to ^ — a known result obtained also from the 
non-abelian integrals for the SU(2) gauge group [9] - [11]. 

In order to study the large cutoff limit more quantitatively we have performed a 
number of extrapolations assuming various asymptotic behaviors of the regularized index 
in n^inax • For example. Fig. 17 shows the asymptotic value extracted with the aid of the 
diagonal Fade approximant P[4:,4]{n% at various T. Two lines correspond to even and 
odd n^max cases, which were independently analyzed . Both extrapolations are stable 
and consistent in the range 1 < T < 5. This result strongly suggests that the infinite 
UBmax limit of thc regularized index is time independent. Moreover, the limiting value 
is nicely consistent with the above 1/4 (also shown in the Figure). Extrapolations with 
power series in different variables lead to similar conclusions. We expect to accumulate 
new data with yet higher cutoffs. This would allow to extend stable extrapolations to 
larger range of T and possibly distinguish between various asymptotic forms tested so far. 

5 Summary and outlook 

The next step in studying a family of supersymmetric Yang-Mills quantum mechanics has 
been completed. Thc above models appear in many areas of theoretical physics, beginning 
with the soluble D = 2 systems, through the small volume, lattice studies of the D = 4 
QCD, and finally ending on the D — 10 models of the M-theory. We are now somewhere 
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in the middle of this hst. 

The new approach presented here leads to the precise study of the rich structure of 
the D — A system, which already has some features of the D = 10 model. With rotational 
invariance taken fully into account, the Hilbert space splits into channels of conserved 
angular momentum and fermionic number. This eliminates the brute force diagonalization 
of large matrices. Second improvement is brought by the generalization of the recursive 
scheme of computing matrix elements while gradually increasing the harmonic oscillator 
basis [23] . Present results fully confirm and extend findings of the first paper where the 
whole program was originated [6] . 

The system has both discrete and continuous spectrum which coexist at the same 
energies. This rather unusual feature was expected for a long time as a consequence of the 
supersymmetric interactions with flat directions. Now however, more precise statements 
can be made. While the discrete, locahzed states exist in all {nF,j) channels, the non- 
localized ones appear only in the central (with respect to the particle-whole reflection) 
supermultiplets and only for even angular momenta. 

The numerical part of the method requires limiting the Hilbert space. We take as 
a cutoff the maximal number of quanta of all bosonic harmonic oscillators, ub max- The 
present approach allows to reach such a large cutoffs that the lower part of the discrete 
spectrum has practically converged to its continuum (i.e., the infinite cutoff) limit. 

On the other hand the eigenenergies from the continuous spectrum literally never 
converge to their continuum values. Instead, they all fall to zero with increasing cutoff. 
In fact, this is precisely the property allowing a clear distinction of the two spectra, 
cf. Figs 2 and 5. The physical energies of the non-localized states arc coded in the 
rate of fall of the above levels with n^max- The particular scaling which governs this 
behavior was discovered some time ago [37] and is well confirmed with present data. It 
is an important tool in extracting any observable related to the non-localized states. In 
particular it allowed us to establish the dispersion relation for the scattering states in the 
{np = 2,j = 0) and {np = 3,j = ^) channels. 

Super symmetry is broken by the cutoff. Again however, with currently available val- 
ues of UBmax, we observc clear restoration of SUSY which manifests itself in many ways 
in the discrete spectrum. First, the energy levels from different channels, related by su- 
pcrsymmetry, coincide to high accuracy, cf. Fig. 6 . Second, our approach allows to form 
and analyze the supersymmetric images of arbitrary eigenstates. This led to the con- 
struction of the rotationally invariant supersymmetry fractions which provided a simple 
identification of SUSY partners. A number of lower supermultiplets was identified for a 
range of angular momenta, see Tables 9 and 10. Interestingly some of the supermultiplets 
are degenerate, see Fig. 13. We do not know a symmetry (if any) responsible for this 
additional degeneracy. The mixing angles are stable with respect to changing the cutoff. 
Their actual values, however, may be an artefact of our regularization. 

A third method to see restoration of supersymmetry is provided by the Witten index. 
It is clearly flattening as a function of euclidean time when we move towards bigger cutoffs 
which are now available. This shows that the cancellations between supersymmetric 
partners becomes more and more efficient, also globally. At infinite cutoff contribution 
from localized states would be exactly zero. Supersymmetric vacuum and other non- 
localized states should give the final non- integer value 1/4 for the gauge group considered 
here. We see now much stronger evidence for this behavior than in the first attempts. 

In the continuous sector of the theory the situation is more difficult and challenging. 
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Although the scaUng expected from the one dimensional free case has been confirmed, it 
should be studied now more extensively, also for higher angular momenta. Identification 
of the supermultiplets is more delicate and remains to be done. Similarly revealing a 
signature of the SUSY vacuum requires further study and yet higher cutoffs. On the 
other hand current precision allows to address more advanced problems like the scattering 
[19, 20]. We are looking forward to work out some of these questions. 

Progress towards higher n^niax for the D = 4 system is limited by computer time. 
The recursive algorithm is presently implemented in Mathematica. We re-implemented 
some sections of the algorithm in C++, obtaining a 100-fold increase in speed; we plan 
to complete the C++ implementation and to improve the present computation. 

Altogether the present approach works rather well. As such it provides one route of 
attacking higher dimensions. The gain from exploiting fully SO(D— 1) invariance and 
restricting ourself to a particular representation of SO(D— 1) should overcome the huge 
sizes of bases in higher dimensions. Generalization to = 10 requires in particular 
construction of the Clebsch-Gordan coefficients for the SO (9) group which is a reasonably 
tedious but a well defined exercise. Some work in this direction has already begun. 
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A Useful identities involving Clebsch-Gordan coeffi- 
cients, 3j symbols, and 6j symbols 

With the usual phase conventions, the Clebsch-Gordan coefficients C^/^-^^ are real and 

C'mim20 = _|_ I " '^"^i+"^2,0; (47) 

the completeness formulae read 

j,m mi,m2 

The Clebsch-Gordan coefficients can be written in terms of the Wigner 3j symbols as 



the 3j symbols enjoy the symmetry properties 



32 Jl J3 j ^ ^_lYi+j2+j3 I Jl J2 J3 \ /gQN 
777-2 nT'l ) V ^1 "^2 7773 / ' 



-Till —772.2 —Tn^ j [mi 7712 ^3 



(51) 
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we also need the formula [40] 

E (-1 



J4+j5+j6-m4-m5-m6 I Jl J5 ] I -^2 -^6 W J4 J5 J3 

\ nil —1715 me ) \ 777.4 —me I \ —7774 7715 777.3 



1714 ,1715 ,me 

^ f ji 32 h \ f Ji 32 h 

I 7711 7712 7773 i I J4 J5 J6 



(52) 



where the term in braces is the Racah 6j symbol. Eqs. (49), (50), and (51) imply the 
"exchange" formulae 



r<32ji3 — ( _-[\ji+j2-j ryjij2j r<hjj2 _ ( 1 \i-j2-m1 1 ^3 2 + 1 ^ ji j2 j (t-o\ 

^m2m\m \ ^) ^m\m2mi ^mimm2 \ -^J V 2j + 1 "^i"*2m' \'-''-'J 

and the "inversion" formula 

ri32 h 3 — ( _'[ \3i+32-3 r<3i32 3 ■ (r:\A\ 

—1712 —mi —m \ ^) ^mi_m2mi \'~'^/ 

Eq. (52) implies 



'.36 

'nrumzmi ^ mAmems ^m5m2m6 
m4 ,m5 ,me 



= (-l)^-^+^3+.4+.5 ^(2j, + l)(2j6 + l) I ^.^ ^.J I C^^^^;:„3. (55) 

B Computation of matrix elements of products of 
operators 

We wish to exploit rotation invariance to reduce the number of matrix elements which 
must be computed and stored; our main tool is the Wigner-Eckhart theorem: 

(j',7n';z'|<Ljj,m;^) -^^^^C^r:,;L'{3'■,^'\\Ol^\\3;^), (56) 

where {j';i'\\Ojl^\\j;i) denotes a reduced matrix element, independent on m', mi, and 777. 

A first formula regarding reduced matrix element can be easily obtained applying Eqs. 
(53) and (54): 

{3';no,M^) = i-iy^^^-'' O,,^ = (-1)^+"^ O,,-^- (57) 

A second formula deals with the product of two operators 0jj^m.i Oj^^,^^ with given 
bosonic and fermionic number: by decomposing the of the product in components with 
definite and then applying Eq. (55), we obtain 

{3'-,i'\\io^]^,o^^)jr,i) 

3+3'+3l+32 JoA I 1 J -^3 3 1 •'||/oW|U-"-.-"\ 



I 31 32 ) 

Applying Eqs. (57) and (58), we never need to deal explicitly with 77i's and Clebsch- 
Gordan coefficients, and the computations are much more efficient. 

(Anti)commutators of operators are dealt with in a very similar way. 
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C Implementation of the algorithm 



We implemented our algorithms in the symbolic manipulation program Mathematica. 

We first compute tables of all needed (anti) commutators and decompositions of prod- 
ucts of two trilinear operators into products of three bilinear operators: we define an 
explicit and univocal representation of a generic operator in terms of a. oJ, /, f'^. in a 
"canonical" order; using this representation, we compute explicitly the desired operators 
and decompose them in the appropriate basis of gauge-invariant operators. The compu- 
tation so far is exact, and the coefficients are square roots of rational numbers. Many 
checks are performed: besides verifying the rotational properties, we check explicitly Eqs. 
(26), (27), (29), (31), and (32). This step requires moderate computer resources. Once 
the tables are computed, the explicit representation of the gauge-invariant operators is 
no longer needed. 

A separate program reads in the tables and implements the orthonormalization and 
recursive computation of scalar products and matrix elements described in Sects. 3.4 and 
3.5. The formulae given in the two sections, together with Hcrmiticity and decomposi- 
tions of products of two trilinear operators, are more than sufficient to reduce any matrix 
element to matrix elements involving a lower number of elementary creation and annihi- 
lation operators; in many instances, more than one reduction is available, and the choice 
can affect performance very strongly. 

For performance reasons, it is crucial to "remember" the values of all matrix elements 
already computed, and to save them periodically into a file to be able to restart the 
computation. Again for performance reasons, we choose to represent matrix elements as 
double-precision floating point numbers rather than as exact algebraic numbers. 



The algorithm described in Sec. 3 generates bases in each channel {nF,j), recursively in 
ub, by applying all operators listed in Sees 3.1, 3.2. Then the Gram-Schmidt orthonor- 
malization selects maximal set of linearly independent states. Their numbers are quoted 
in Tables 12-15. 

Equivalently, dimensions of the above-mentioned subspaces can be derived by classi- 
fying all independent tensor structures contributing to each channel at given ub- This 
provides an additional check of our program and prepares the ground for subsequent gen- 
eralization to higher dimensions and higher gauge groups. Here are few examples for 
various np- 

D.l AU riF = states. 

For even ub, every gauge- invariant state can be obtained by applying to the vacuum a 
combination of the gauge-invariant creation operators 



Since there are 6 independent Aj^, and states created by different products of A]/^, apart 
from permutations, are linearly independent, the total size of the basis with even number 
of bosons is 



D Sizes of bases 



(59) 
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D.2 np = 0,j = 0,2. 

All gauge invariant and spherically symmetric states can be obtained by combining the 
traces of the powers of the basic gauge invariant bihnear creator (59). Since A is a three 
by three matrix, its Cayley-Hamilton equation is third order, hence only traces of first 
three powers of A are independent. It follows that the number of independent states with 
ub quanta equals to the number of monomials of the nB/2 order which can be made from 
Tt A, Tr(yl^) and Tr{A^). Therefore it is given by the number of partitions P(ns/2|1, 2, 3) 



of ns/2 into elements smaller than 4. 

States with odd ub can be generated by acting with the only odd creator A^, Eq. (24), 
on the even basis. Therefore 



This explains the even-odd regularities in the first column of Table 12. Since IS a 
scalar Eq.(61) holds for arbitrary angular momentum j and consequently also for global 
number of states, cf. Table 2. 



J\f{0, 0, Ub) = P(nB/2|l, 2, 3), - even, 



(60) 



J\f{j = 0, riF = 0, Ub) = A/'(0, 0, ub - 3), ub - odd. 



(61) 
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riB 
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1 
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1 
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1 





1 
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1 







































1 


8 


4 


2 





2 
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21 


29 


5 


1 





1 



































6 


35 


6 


3 





3 


1 


2 





1 
























56 


91 


7 


2 





2 





1 































21 


112 


8 


4 





5 


1 


4 


1 


2 





1 




















126 


238 


9 


3 





3 


1 


2 





1 



























56 


294 


10 


5 





7 


2 


6 


2 


4 


1 


2 





1 
















252 


546 


11 


4 





5 


1 


4 


1 


2 





1 























126 


672 


12 


7 





9 


3 


9 


3 


7 


2 


4 


1 


2 





1 












462 


1134 


13 


5 





7 


2 


6 


2 


4 


1 


2 





1 



















252 


1386 


14 


8 





12 


4 


12 


5 


10 


4 


7 


2 


4 


1 


2 





1 








792 


2178 


15 


7 





9 


3 


9 


3 


7 


2 


4 


1 


2 





1 















462 


2640 


16 


10 





15 


5 


16 


7 


14 


6 


11 


4 


7 


2 


4 


1 


2 





1 




1287 


3927 


17 


8 





12 


4 


12 


5 


10 


4 


7 


2 


4 


1 


2 





1 











792 


4719 


18 


12 





18 


7 


20 


9 


19 


9 


15 


7 


11 


4 


7 


2 


4 


1 


2 





1 2002 


6721 



Table 12: Number of S0(3) multiplets with np = and fixed j and Ub- Ns is the number of basis vectors with given number of 
bosonic quanta, hb, while E gives the cumulative size up to ub- 



The j = 2 states can be generated from the empty states by replacing one of the two 
traces Tr(74) or Tr^A^) by the symmetric traccless tensor formed from A or A^. Note 
that such a tensor formed from is already dependent on the lower powers of A. This 
is again the consequence of the Cayley-Hamilton equation: only the trace of the is 
independent since it is equivalent to one of the coefficients of the C-H equation (namely 
to the determinant). One can therefore count the j=2 states as follows: for each even ub 
take all partitions contributing to Eq. (60), replace in each partition one element, e.g., 2 
by its indexed counterpart 2*^^. This produces a monomial with Tr(A^) — * (A^)**^ — A'^5^'^ 
which generates one j — 2 state. Repeat this procedure for all different elements in a 
partition omitting value 3. Total number of states equals to the number of such indexed 
monomials. This procedure indeed reproduces sizes listed in the third column of Tablel2. 
Yet simpler counting can be formulated recursively: states with ub quanta can be obtained 
by acting with A^'' on the j — 0,nb — 2 basis and independently by acting with {A^y'' on 
the j — 0,nB — 4: basis. This gives the recursion relation 

^^{Q, 2, Ub) = 7\A(0, 0, ns - 2) + M{Q, 0, - 4), (62) 

which explains the j — 2 column of Tablel2. 
D.3 np = l,j = 1/2,3/2. 

The lowest gauge invariant state in this sector must contain one boson and is created by 
fl'^al^ = {faY^ from the empty state. States with j = 1/2 and j = 2/3 are generated by 
suitable projections (/a)i/2 and {f 0)3/2, where (...)j means summing over a and i indices 
with appropriate Clebsch-Gordan coefficients. For odd ub one then combines powers 
of traces of A.A"^ and A^ from previous case with three independent^ fermionic creators 
(/'2)i/2i ifo.A)i/2, (/a.A^)i/2 to get all states with j = 1/2 ^. One can generate all states 
of the hb basis recursively by acting with {f 0)1/2 on the np = 0, j = 0, — 1 basis; with 
{fo.A)i/2 on the np = 0, j = 0,nB — 3 basis; and {fo.A^)i/2 on the np = 0, j = 0,nB — 5 
basis. This imphes the relation 

A/'(l, 1/2, ns) = ^^{Q, 0, - 1) + A/'(0, 0, ns - 3) + X{Q, 0, - 5), (63) 

which gives the first column of Table 13 in terms of Table 12. 



^ Again Cayley-Hamilton equation for A limits a number of independent creators with rii? = 1. 
^A "." denotes contraction of an adjoint S0(3) indices. 
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4 
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5 


4 


2 
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126 


210 


6 


4 


6 


5 


4 


2 
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126 


336 
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10 


10 
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2 
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336 


672 


8 


6 


10 


10 


9 


6 


4 


2 


1 

























336 


1008 


9 


9 


15 


16 


16 


13 


10 


6 


4 


2 


1 




















756 


1764 


10 


9 


15 


16 


16 


13 


10 


6 


4 


2 


1 





















756 


2520 


11 


12 


21 


24 


25 


22 


19 


14 


10 


6 


4 


2 


1 
















1512 


4032 


12 


12 


21 


24 


25 


22 


19 


14 


10 


6 


4 


2 


1 

















1512 


5544 


13 


16 


28 


33 


36 


34 


31 


25 


20 


14 


10 


6 


4 


2 


1 












2772 


8316 


14 


16 


28 


33 


36 


34 


31 


25 


20 


14 


10 


6 


4 


2 


1 













2772 


11088 


15 


20 


36 


44 


49 


48 


46 


40 


34 


26 


20 


14 


10 


6 


4 


2 


1 








4752 


15840 


16 


20 


36 


44 


49 


48 


46 


40 


34 


26 


20 


14 


10 


6 


4 


2 


1 









4752 


20592 


17 


25 


45 


56 


64 


65 


64 


58 


52 


43 


35 


26 


20 


14 


10 


6 


4 


2 


1 




7722 


28314 


18 


25 


45 


56 


64 


65 


64 


58 


52 


43 


35 


26 


20 


14 


10 


6 


4 


2 


1 





7722 


36036 



Table 13: Number of S0(3) multiplets with np = 1 and fixed j and Ub- Ns is the number of basis vectors with given number of 
bosonic quanta, hb, while E gives the cumulative size up to ub- 



For j = 3/2 new fermionic creators can be constructed beginning with = 3. In this 
case there are two old creators (/a)3/2 Tr(74) and {fa.A)3/2, corresponding to couphng 
3/2 ® — > 3/2, and (1/2 (g) 1) 2 — > 3/2 ®. However there are two ways to reahze the last 
coupling. Hence there must exist an independent creator (fa. A) 3/2, where A^"' denotes 
gcnerically the bilinear form of bosonic creators convohitcd with all appropriate Clebsch- 
Gordan coefficients. Now, to generate all independent states with np = l,j = 3/2 and 
Ub bosons we have to act with (/a)3/2 on the np = O.j = 0,72^ — 1 basis, with {fa.A)3/2 
and {fa.A)3/2 on the np = 0,j = 0,nB — 3 basis, with (/a.A^)3/2 and {fa.A.A)3/2 on the 
np = 0, j = Q,nB — 5 basis, and finally with {fa.A.A.A)3/2 on the np = 0, j = 0,nB — 7 
basis. Therefore 

1/2, ns) = 7\A(0, 0, nB-l)+2(7\A(0, 0, ns-3)+A/'(0, 0, ns-5))+A/'(0, 0, ns-7). (64) 

This explains the second column of Table 13. 
D.4 np = 2,j = 0. 

A general two-fermion creator can have three covariant forms: (a) an SO (3) scalar F ~ 

fb^^'^fb ' (^) S0(3) scalar symmetric in color indices '^F^c = fl^^"^ fc~^^^—fb ^^^/c^^^j 
and (c) an SO (3) vector antisymmetric in color indices ^F^^ corresponding to a cou- 
pling: 1/2(8)1/2— >1. To construct gauge invariant creators from (2) define the matrix 
{^FAy'' Fbcttl^al''. Independent symmetric creators are 

Tr(^FA), Tr(^FAA), Tr{^ FA.A^) and F. (65) 

Antisymmetric creators are constructed from the matrix ("^^^4)*' =^ F^^a]^ a]}' e^^K Cayley- 
Hamilton equation allows to construct two independent creators in this case 

Tr(^FA), Ti{^FA.A). (66) 

Then the np = 2,j = 0,basis with quanta can be obtained by acting with all these 
creators on the appropriate j = bases with lower ^. As a consequence 

J\f(2, 0, ns) = J^(0, 0, ns) + 2(X(0, 0, ns - 2) +X(0, 0, ns - 4)) +X(0, 0, 71^-6). (67) 

Which is readily satisfied by the first column of Table 14. 
D.5 np = 3,j = 1/2. 

It is left as an exercise for the reader to prove why the first column of Table 15 satisfies 
the following recursion for odd ns 

J\f{3, 1/2, ne) = 7^(0, 0, - 1) + 2A^(0, 0, tib - 3) + (0, 0, tib - 5) + 3J\f{0, 0, ns - 7). 

(68) 



^The trace in A is linearly dependent with the the first creator and there is no gauge invariant, j=l, 
combination of a's 

^This decomposition is analogous to that of van Baal [35] . 
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15 
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39 


32 


34 
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18 
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4025 


10 
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21 


49 


36 
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27 


27 


12 
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6146 


11 


21 


30 


57 


50 


57 


42 


38 


23 


18 
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6 


2 
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3234 


9380 


12 


28 


30 


69 


55 


71 


49 


51 


29 


27 


12 


11 
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3 
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13566 


13 


28 


40 


79 


72 


86 


68 


67 


46 


39 


23 


18 
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19614 


14 


36 


40 


93 


78 


103 


77 


84 


55 


53 


29 


27 


12 


11 
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3 
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27210 


15 


36 


52 


104 


98 


121 


101 


104 


78 


71 


47 


39 


23 


18 


9 
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2 
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10530 


37740 


16 


45 


52 


120 


105 


141 


112 


125 


90 


90 


57 


53 
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27 


12 


11 


3 


3 


12915 


50655 



Table 14: Number of S0(3) multiplets with np = 2 and fixed j and ub- Ng is tfie number of basis vectors with given number of 
bosonic quanta, hb, while E gives the cumulative size up to ub- 
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11 
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39 
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42 
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11 
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12 


39 
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58 


44 


29 
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11 
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13 


42 


77 


93 


100 


93 


83 


66 


50 


34 


23 


13 
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14 


52 


96 


119 


131 


127 


118 


100 


81 


60 


44 


29 


19 


11 
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2 
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11232 


36600 


15 


54 


102 


126 


139 


137 


128 


109 


91 


68 


50 


34 


23 


13 


8 


3 


1 







12480 


49080 


16 


66 


123 


156 


176 


177 


171 


153 


132 


106 


83 


60 


44 


29 


19 


11 


6 


2 


1 


18900 


67980 



Table 15: Number of S0(3) multiplets with np = ?> and fixed j and Hb- Ng is tlie number of basis vectors in all angular momentum 
channels with given number of bosonic quanta, n^, while E gives the cumulative size up to n^. 
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